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Abstract 

A self-consistent and universal description of friction and diffusion for Brow- 
nian particles (grains) in different systems, as a gas with Boltzmann colli- 
sions, dusty plasma with ion absorption by grains, and for active particles 
(e.g., cells in biological systems) is suggested on the basis of the appropriate 
Fokker-Planck equation. Restrictions for application of the Fokker-Planck 
equation to the problem of velocity-dependent friction and diffusion coeffi- 
cients are found. General description for this coefficient is formulated on the 
basis of master equation. Relation of the diffusion coefficient in the coordinate 
and velocity spaces is found for active (capable to transfer momentum to the 
ambient media) and passive particles in the framework of the Fokker-Planck 
equation. 

The problem of anomalous space diffusion is formulated on the basis of 
the appropriate probability transition (PT) function. The method of partial 
differentiation is avoided to construct the correct probability distributions for 
arbitrary distances, what is important for applications to different stochastic 
problems. Generale equation for time-dependent PT function is formulated 
and discussed. 

Generalized friction in the velocity space is determined and applied to de- 
scribe the friction force itself as well as the drag force in the case of a non-zero 
driven ion velocity in plasmas. The negative friction due to ion scattering on 
grains exists and can be realized for the appropriate experimental conditions. 

PACS: 52.27.Lw, 52.20.Hv, 05.40.-a, 05.40.Fb 
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I. INTRODUCTION 



Interest in Brownian dynamics is conditioned by a large variety of applications: granular 
systems, including dusty plasmas, various objects in biological systems, physical-chemical 
systems, et cetera. For open systems and for systems with non-elastic processes the rela- 
tions between the friction and diffusion coefficients and even the correct specific forms of 
the Fokker-Planck equation are not still established. Recently the friction and diffusion 
coefficients as functions of the grain velocity V were derived for dusty plasmas [1]. Due to 
ion absorption by the grains, the friction coefficient can become negative [2] for the sim- 
plest model of dusty plasmas with ion absorption without atom regeneration on the grain's 
surface. This implies that the problem of negative friction has to be considered for more 
complicated models and realistic situations in dusty plasmas, in particular taking into ac- 
count the grain's mass conservation for, which is essential for large times of the process, 
in spite of big difference between ion and grain masses. In the present paper we consider 
the specific forms of the probability transition (PT) [3] for Boltzmann-type and and non- 
Boltzmann-type of collisions to calculate and compare the velocity dependent friction and 
diffusion coefficients in different systems, including gases, dusty plasmas and active particles. 
Our main goal in this paper is determination of the velocity-dependent friction and diffusion 
coefficients for arbitrary grain velocity and consideration of diffusion in the coordinate and 
velocity spaces for different systems and conditions. 

In Section 2 the Fokker-Planck equation and self-consistent description of velocity- 
dependent friction and diffusion are discussed in details on basis of the PT function (PT). 
In Section 3 we consider in detail the probability transition function for Boltzmann-type 
collisions, in particular, for the Boltzmann's spheres. Friction and diffusion coefficients are 
calculated on the basis of the Fokker-Planck equation. The restriction for apphcation of the 
Fokker-Planck equation to the fast motion of grains for the Boltzmann's type of collisions 
and necessity to use the appropriate master equation with prescribed PT function is shown. 
The general formulae for the velocity-dependent friction and diffusion coefficients in the 
velocity space arc established on the basis of master equation. In Section 4 the absorption 
collisions for ions by grains is done for dusty plasmas with a high grain charge. In the frame- 
work of the simplest model with a fixed mass of grains the phenomena of negative friction 
is found. The friction coefficient changes sign for some value of the grain velocity. Mani- 
festation of negative friction due to ion absorption in dusty plasmas is the consequence of 
the model. Opportunity to realize this process in the experiment can be considered on more 
elaborated model with grain mass conservation. In Section 5 we consider diffusion in the 
coordinate space on the basis of the Fokker-Planck kinetic equation. The general relation 
between the diffusion coefficient in the coordinate space and velocity-dependent diffusion 
and friction coefficients in the velocity space are found. Relation of these results with the 
known limiting cases is established. The results are applied to the self-motion of grains (e.g. 
cells) in biological systems. In Section 6 the different types of anomalous diffusion in the 
coordinate space are considered on the basis of the appropriate master equation. As partic- 
ular case we describes the Levy-flights, which are very important for many applications. We 
also formulated the generalized diffusion equation for the time-dependent PT function in the 
coordinate space, investigated it for some particular cases and suggest to apply for subdiffu- 
sion processes. In Section 7 we shortly describe the generalized friction, which appears when 
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some other vectors, except grain velocity, are essential to determination of PT function in 
the velocity space. The application of such consideration is done in Section 8, where the 
ion flow with some driven velocity (which is the particular case of such additional vector 
in PT function) scatters on grains. Due to this process the total friction can be negative 
for small grain velocity and provides the mechanism for acceleration and heating grains in 
dusty plasmas. The conclusions are presented in the Section 9. 



II. FOKKER-PLANCK EQUATION AND SELF CONSISTENT 
VELOCITY-DEPENDENT FRICTION AND DIFFUSION 

The friction and diffusion coefficients P{V), {V) and D±{V) on the basis of the Fokker- 
Planck equation are determined via PT function w(P,q) [3] as 

W = / rf^q(Pq)^P,q), (1) 

= (1/2M^P2) J d^q(Pq)^^(P,q). (2) 

and 

D±iV) = _ l^^,p, J rf^q [pV - (Pq)1 ^(P, q)- (3) 

Here (Pq) is the scalar product in velocity space of the dimension s and P = MV, where 
M and V are respectively the mass and velocity of the grain. 

We determine the PT function w{P,q) = w(P,P' — P) (q = P' — P represents the 
momentum transfer for the considered process, e.g. collision) in the same way (convenient 
for consideration of the Fokker-Planck equation) as was done in [5] : 

Ig{P,t) = y"dq{w(P + q,q)/,(P + q,t) 

-^P,q)/,(P,i)}, (4) 



dt 



The corresponding to Eq. (4) Fokker-Planck equation can be written as 

df,{P,t) _ d 



dt dPi 



P,/3(P)/,(P) + A (A,(P)/,(P)) 



(5) 



The PT function w(P,q) is directly connected with the probability transition u'o(P,P') 
in the Chapman-Kolmogorov equation (master equation, [4]) wo(P',P) = w(P,P — P') = 
w(P,q): 



dt 



= j dP'{«;o(P,P')/.(P',t)- u.'o(P',P)/,(P,t)}. (6) 



It follows from Eqs. (1), (2) that the friction and diffusion coefficients always have to be 
calculated self- consistently. For the Boltzmann-type collisions between the light particles 
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(conventionally atoms) and grains with masses m and M (m ^ M), respectively, the PT 
function ^^(P, q) and the coefficients f3{V) and D(y) can be found for arbitrary cross-section 
(as will be done in Section 2). In Section 2 we also investigate the applicability of the Einstein 
relation to the Boltzmann-type collisions on the basis of the Fokker-Planck equation with 
velocity-dependent coefficients. 

For some other type of collisions the problem of a single-valued form of w(P, q) exists, 
because the exact kinetic equation for such collisions is not known in many cases. The 
Einstein relation between friction and diffusion of grains can be easily violated for such type 
of non-Boltzmann collisions, even for very low grain velocity. 

The solution of the stationary Fokker-Planck equation for the isotropic distribution func- 
tion of the particles can be written [3] as: 



c 



cxp 



V 



— \ dvv 



D\\{v) 



(7) 



where C is a constant, providing normalization and D\\ is the longitudinal part of the diffusion 
tensor: 



A,(P) = ^D\\{P) + (% - ^) D^{P). 

The function (3*{V) is related with I3{V) by the equality: 

s- 1 



(3*{V) = (5{V) + (L'ii(y) - D^{V)) . 



From Eq. (7) follows: 



1 dL 



V dV Vfg dV 



(8) 



(9) 



(10) 



Here the function fg is the stationary grain distribution function. 



III. BOLTZMANN-TYPE COLLISIONS: LIMITATIONS OF THE 
FOKKER-PLANCK EQUATION. MASTER EQUATION APPROACH 

For the Boltzmann-type collisions the stationary distribution fg{V) is a Maxwellian equi- 
librium fgo{V) and Eq. (10) has the form 

If D is a function of V the Einstein relation (5 = MD/T is violated if the term with derivative 
in Eq. (11) is not small with respect to the parameter m/M. Only in that special case 
deviation from the Einstein relation in Eq. (11) can be neglected, according to the general 
conditions of applicability for the Fokker-Planck equation. For Boltzmann collisions it is 
possible to show that the PT function ^^(P, q) (with P = Pa = MV) can be represented 
as: 



4 



WB{P,q) 



(iu 5 ( uq — 



2/1 



da 
do 



arccos (1 — 



),u 



A(u + V-vo) 



(12) 



Here ^ is the reduced mass of the particles of species a and b (with masses M and m, 
respectively), da /do is the differential cross-section for Boltzmann-type scattering and is 
the velocity distribution function for the light particles, which we assume to be Maxwellian 
with, for generality, some shifted velocity. It means small particles can have some driven 
velocity vq. 

For the case of hard spheres, when da /do — a^/A and the distribution function of the 
small particles fb — fbo is Maxwellian , Eq. (12) can be written as: 



<(P,q) 



Aqjj? \2nTb 



1/2 



cxp 



rrih 
'2T 



Vq 

Q 



2/i 



(13) 



where a is the reduced radius of the colliding spheres. Let us calculate the friction and 
diffusion coefficients on the basis of general equations (1),(2) and the PT function (13) for 
hard spheres. It is easy to show that the friction coefficient can be represented via the 
function of parabolic cylinder D_3 in the form: 



+ 



uP'x 



D_ 



(14) 



where uj = V/vt and vt — yjT/m is the thermal velocity of the small particles. The value 
of the friction coefficient for the grain in the hard spheres media for y = is equal: 

(15) 



/3(0) = 8V27r/ina^UT/3M. 
Likewise for the friction coefficient we derive 



D(V) - D(0) J , J=^J dxx^ exp 



+ 



2 

UJ X 



D_4{—U!X). 



2 ' 4 

For V = the coefficients D{0) and /5(0) are related (/x ~ m) by Einstein's law: 

D(0) = iiT(3(0)/mM. 



(16) 



(17) 



For the opposite case of very high velocity V ^ vt the formal expressions for the coefficients 
P and D, based on self-consistent relations (1), (2) and the PT function for Boltzmann-type 
collisions (12) are: 



P{V » vt) = ^a'nV, 



(18) 



DiV » vt) = P{V » .^)^ = ^^a'nV^ (19) 

The similar expressions have been found in [3] (although here we corrected the numerical 
multipliers order of unit) and evidently violate the Einstein relation. As it follows from 
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Eq. (13) the conditions of applicability of the Fokker-Planck equation cannot be restricted 
only by the usual inequality (e.g. [5]) q ^ P. To find the full limitations for transition from 
the Boltzmann equation to the Fokker-Planck equation, let us consider the conditions for 
expansion of the function w(P + q) by q in Eq. (4) for the simple particular case described 
by Eq. (13): 



(P + q, q) = «;^(P, q) exp 



rrib /Vq 
'2T [ q 



_q_\ 2q 
2fiJ M 



mb q 



(20) 



To obtain the Fokker-Planck equation from the master equation (4) we have to expand 
the exponent in Eq. (20). Therefore besides the general condition q <^ P the inequalities 
I V |<S v^M/q and q'^/TM <S 1 have to be fulfilled simultaneously. An even more stronger 
limitation for applicability of Fokker-Planck equation for Boltzmann- type collisions | V | <^ 
T/q follows from the fact that the stationary state of the distribution function for grains 
(without external fields) is Maxwellian. This means that for any finite q there is some 
limitation on the values | V |, which can be described by the Fokker-Planck equation. In 
particular, according these inequalities the limit T — > cannot be realized for the Boltzmann- 
type collisions on the basis of the Fokker-Planck approximation simultaneously with finite 
I V I -values. In other words, the expressions for the friction and diffusion coefficients 
(18), (19) for a Boltzmann gas for the velocities, which are not satisfied the inequalities, 
mentioned above, are in contradiction with Eq. (11). 

Because the Fokker-Planck equation itself is not valid in Boltzmann gas for high ve- 
locities, Eqs. (5), (11) is not applicable and our aim is to find the friction and diffusion 
coefficients for arbitrary velocities on the basis of the Chapman- Kolmogorov equation (6). 
For this equation in the absence of the external fields the stationary distribution for grains 
is, as easy to see, the Maxwellian if the transfer probability function is calculated with the 
Maxwellian distribution of the light particles (13). Equation (6) permits to determine in a 
physically justified way the friction and diffusion coefficients for arbitrary particle velocities 
in the same form (1), (2) as it was found from the Fokker-Planck equation, in spite the 
equation itself is not valid anymore. 

In order to show this let us integrate equation (11) with the multiplayer and with the 
product PaPp- Then we arrive to the equations for the momentum: 



rf(P) 
dt 



- j dP/,(P) j c^P'«;o(P',P)(P-P') = 
j dPf,{P) J dq«;(P,q)q= - J dPf,{P)Pf3{P). 



(21) 



Here (P) is the momentum, averaged by the distribution function fg and /9(P) is exactly 
the same as in Eq. (1). Equation for the averaged product (P^P^) is: 



dt 



- y dPfg{p) J dP' (p,p^-p',p'^)wo(p',p)= 



(m2D||(P) - P'PiP)) + 2M^ - 



PgPp 
p2 



D±{P) 



2M^ J dPfg{P) 



Dq,{p) - 



(22) 
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The function -D||(-P) D±{P) are the same as in Eqs. (2). (3). which were found from the 
Fokker-Planck equation, ahhough this equation is strictly not apphcable and we derived 
Eqs. (21), (22) on the basis of master equation (6). 

If we recalculate the time derivatives for the averaged momenta on the basis of Fokker- 
Planck equation (5) we would find the same form of the resulting equations (21,22). Natu- 
rally, the function fg in both cases is different and satisfies to the relevant (Chapman-Enskog 
or Fokker-Planck) equations. It means that determination of the friction and diffusion co- 
efficients by Eqs. (l)-(3) is general for the system of a two components, if the distribution 
function for the light component is prescribed. In this way we express the equation for the 
momentum of distribution function with the general velocity-dependent friction and diffu- 
sion in the velocity space. Of cause, the relation (10) is not valid anymore, nor are the 
Einstein relations for the velocity dependent coefficients. This consideration justifies and 
generalizes the statements of [3]. Interpolation relation between the friction and diffusion 
coefficients, which is vahd for all velocities for the considered case can be written in the 
form: 




which corresponds with the Einstein relation only for the case that the velocities are much 
smaller than thermal one V <^ vt- In this sense a discussion of the relations between the 
velocity-dependent friction and diffusion coefficients on the basis of the Fokker-Planck equa- 
tion [6] is relevant for non-Boltzmann collisions; for Boltzmann-type collisions this discussion 
seems artificial, if there is not some specific additional parameter, which makes the condition 
q <^ P sufficient for using of the Fokker-Planck equation, as it is for the case of collisions in 
a plasma with Coulomb interaction. 

At the same time in parallel with the interpolation relation (23), which, naturally, is 
not single- valued, the exact integral relation between the velocity dependent friction and 
diffusion exists as it follows from Eq. (22): 

J dPfgo{P)D^^{P) = J dPfgoPaP^PiP). (24) 

This relation is the integral generalization of the Einstein's law and transfers to it for the 
case of velocity-independent friction and diffusion coefficients. 



IV. ABSORPTION COLLISIONS 

For absorption collisions (typical, e.g., for dusty plasmas) the momentum transferred to 
the grain is equal to the momentum of the colliding atom (ion in the case of dusty plasmas) 
p. The PT function w{P, q) in this case can be found [3] on the basis of the kinetic equation 
for absorption [7] 



w{P,q)^ fni-flWbs 



P q 

h — 

M m 



P q 

h — 

M m 



(25) 
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For the Maxwellian distribution of atoms and purely geometrical absorption aabs = ^ra^ 
(the particular case of the absorption cross-section for the model of dusty plasmas with the 
average charge of grain Q — > 0), we find for arbitrary V: 

D(V)^2M~^T(3(V). (26) 

For V — this result was found in [1]. Equation (26) is different from the Einstein relation 
because the absorption collision integral is different from the Boltzmann-type. 

For the real dusty plasmas, where the average charge of grain Q is very high (order of 
of the electron charge) the absorption cross-section can be taken in so called orbital 
motion limited (OML). This approach based on the assumption that there is no potential 
barrier for the ions moving towards the grain and the conservation of angular momentum 
and energy is enough to obtain the cross-section aabs = ctoml of absorption for the charge 
Ca by grain with the radius a: 



aoML{Q,Va) = vraM 1 ^ ^1 ^ . (27) 

Here V = |v„ — V| is the relative velocity of the colliding particles. As it follows from 
Eqs. (26), (27) the Fokker-Planck equation is always apphcable to the absorption coUisions 
under consideration if P » g due to smooth dependence of the PT function from the 
momentum P. Calculation of the velocity-dependent friction and diffusion coefficients on 
the basis of PT function (25) and Eqs. (l)-(3) leads to the conclusion, that friction coefficient 
in the model of dusty plasmas with dominant absorption collisions always negative in some 
velocity domain if the parameter of ion-grain interaction P = Qci/aTi [2]. For the range of 
parameters rj = V^/v^^ <^ 1 and 1 < P < 2 the absorption friction coefficient for ions Pi has 
a form: 

A(r;,P)=2A[l-P-^(l-3P)l, (28) 
L 5 J 

and negative for rj < 5(P — 1)/(3P — 1). The coefficient A is equal: 

A = -V27r ( — ) a^riiVTi. (29) 
3 \mgj 



For T) ^ 1 and arbitrary values of P the equation (1) gives: 

1 + 2P 



2r) 



(30) 



It means that for P ^ 1 and large r] the equation /3j(?7iP) = has a root 77i(P) ^ (1-I-2P) /2 ~ 
P. Naturally, for all P > 1 there exists appropriate r){r), which is the root of the equation 

Necessary to mention that it does not mean that the total friction in dusty plasmas can 
be negative, because there are other mechanisms of positive friction (e.e. atomic and ion 
scattering by grains) which provide usually totally positive friction for grains. Manifestation 
of negative friction due to absorption requires the very special condition. In particular we 
have stress that the friction coefficient which we discussed is finite for all values of 
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velocities. Below we consider another mechanism of negative friction, for which always exist 
the experimental conditions. 

We have to stress that Eq. (25) describes the absorption collisions without mass conser- 
vation, otherwise mass of the grains will increase and stationary solution does not exist [8]- 
[10]. In real dusty plasma both cases (stationary and non-stationary regimes) are realized, 
depending from the types of plasma, discharges and parameters. For the case of stationary 
solutions the mass conservation provides by surface recombination of ions and remove of 
the creating atoms back to plasma. Naturally this more complicated situation of absorption 
with further atom emission by grain cannot be described by PT function (25) and requires 
a special consideration. Therefore the mechanism of negative friction considered above has 
mostly illustrative interest as an example of microscopical model for negative friction. 

V. DIFFUSION IN THE COORDINATE SPACE: SELF-CONSISTENT 
DESCRIPTION FOR THE FOKKER-PLANCK APPROXIMATION 

Let us consider the relation of diffusion in coordinate space with diffusion in the velocity 
space on the basis of the Fokker-Planck equation. It means as it was shown above, that for 

the case of Boltzmann gas we have consider a low grain velocity (with the velocity inde- 
pendent friction and diffusion coefficient, when the Fokker-Planck equation is applicable). 
At the same time for the non-Boltzmann type of collisions, such a type as the absorption 
collisions in plasmas or active motion of biological objects (see [3] and the analysis below) 
we can consider the coefficients, which are essentially dependable from the velocity, because 
the Fokker-Planck equation in that case based on the processes of non-Boltzmann type. 
The problem of relation between the diffusion in the coordinate and velocity spaces date 
to Chandrasckhar [11] has been considered in a few papers (e.g. [12]- [14]) for the Fokker- 
Planck equation, but still is not solved self-consistcntly. To find the diffusion Dr in the 
coordinate space we have consider a slightly inhomogeneous in space stationary solution for 
the Fokker-Planck equation (see inhomogeneous variant of Eq. (5): 



^a/,(P,r,t)_ d 



dr dPi 



mP)f,{P,r,t) + A (A,(P)/,(P,r,i)) 



(31) 



For the distribution function one can use the expansion fg{p, r, t) = fg{p, r, t) + /J(p, r, t), 
where /° is has the same shape as the stationary homogeneous solution, but with coordinate- 
dependent density ng{r) and temperature Tg{r). The equation for perturbation fg can be 
written in the form: 

W)^(P,r,t) + (A,(P)^^(P,r,t)) = ^*(^,r) ^ *,(P,r), (32) 
where the functions P and $ are equal: 

1 /"^ 

^^^''^ = "M / dP'P'fSiP'^^)- (34) 
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Because for the slight inhomogeneity one can represent the function ^ as the hnear combi- 
nation of the gradients: 

*(P, r) = A(P) Vn(r) + 7(P)VT(r), (35) 

it is natural to find a solution /J (P, r, t) of Eq. (32) as 

f^(P, r) = //i(P)P . Vn(r) + a^2(P)P • VT(r). (36) 

Taking into account that fg{P) is always proportional to the density of grains n(r) we easily 
determine the functions A(P) and after that /ii(P): 

1 /""^ 



M2P 







"ciP'^exp 



^11(^0 



y dP'V(p" 



(38) 



The function t(P) under the integral in the exponent of Eq. (38) is equal: 

r(P) = 

For self-consistent description according to Eq. (7) we can rewrite Hi{P) as: 



(39) 



A(P') 



M2P Jo i^|,(P')/°(P') 



M3nP 



rf -1 /-oo 

1 ''^'iMn/p. 



The same way can be found /X2(P). 

Because the diffusion current is determined by the relation: 



b = y d^YfliP, r) = -DrV ■ n(r), 



(41) 



we find straightforward in the framework of the Fokker-Planck equation the general relation: 



D 



1 



R 



sM 



d'PP^fii{P) 



1 



d'PPfAP) / dP' 



M^sn.1 —^^s' ' D\\{P')P^{P')Jp, 



dPPfJP), (42) 



where s = 3 for tree-dimensional case and s = 2 for two-dimensional case. The values 

f3{P), -D|| and D± arc determined sclf-consistently by the relations (1),(2) and (3) via the 
probability transition function. Introducing the characteristic momentum Pq, velocity Vq, 
and constant friction Pc and diffusion coefficients by formulae Pq = MVq, /3(P) = PcK-P) 
and D{P) = Dcd{P) let us rewrite Eq.(42) in undimensional form: 



D 



R 



-jr / d^mflir,) / d^, 



^i||(7)/°(7) A 
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duuf^Auj), 



(43) 



where the distribution function fg{r)) (r) = P/Pq) is normahzed to unit: 



'd\\ iv) 



exp 



-A PdCC^ 
Jo ^^rf||(C)J 



1. 



(44) 



The parameter A = (ScVq / D^- It means in general the relation of the diffusion coefficients 
in the coordinate and velocity space (42) can be written as: 



D 



R 



sD, 



(45) 



where we determine the integrals in Eq. (42) as F {(3cVq/Dc). 

Let us consider the particular cases, 
a) For the constant values of (3 and D — D\\ — D±, when the velocity dependence is absent 
we find: 



P = P,T 



1 



M^DP 



dP'\{P') exp 



/3 



2M2D 



(46) 



Let us introduce the effective temperature T* — MD/p. If the stationary distribution 
is Maxwellian = with some temperature T the function A and /ii from Eq. (46) 
transforms to 



-, l^l{P) 



n 



(47) 



where $(a;) is the error function and l/T^ff = (1/T — 1/T*). 

If the Einstein relation MD = T(3 is fulfilled then T* = T and T^ff oo. For that case 
we have Hi — —T/Jq/uM'^D and the diffusion coefficient in the coordinate space related 
with the diffusion in the velocity space in both 2D and 3D cases by the known equality: 



D 



R 



(48) 



For more complicated physical problems, when the Einstein relation is not applicable, as 
e.g. for active particles, necessary to use the general relations derived above. In general 
necessary also to use the self-consistent expressions for the friction and diffusion coefficients 
in velocity space to find correctly the diffusion in the coordinate space, because the last one 
is the functional of both of them, as it follows from the general expressions (42) or (43). 

However to find /? and D self-consistently and determine on that basis the stationary 
distribution function necessary to find the probability transition function, which is known 
on the microscopical basis only for a restricted type of the momentum transfer processes [3] 
(Boltzmann-type collisions, some models of absorption, in particular in dusty plasmas, and 
some types of active motion). For many other cases, where the well developed microscopical 
model of the physical process is absent, the approximate non-self-consistent consideration, 
based on physical intuition, is only one way to obtain the results. 

b) Let us apply now the general relations, found above, to the model of active parti- 
cles (e.g. cells) capable to self-motion. The problem of active friction have been discussed 
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intensively for the last decade on basis of phenomenological approach (e.g. [15]). The micro- 
scopical model of active friction for biological objects was recently suggested in [3]. Gener- 
alization of this microscopical consideration for more complicated and realistic cases, when 
the "driver" vector of the object, introduced in this paper, can rotate and orients on some 
external gradients of the ambient medium, will be considered separately. In that case the 
amplified force, which provides a negative friction is directed to "driver", although passive 
friction is antiparallel to velocity vector, and the active particle can change the direction of 
motion. 

The PT function for the active particles has been found in [3] by considering physical 
processes for momentum transfer from active particles to the ambient medium by internal 
energy loss of the inner energy of the grain (cell). Our main finding is a "singularity" 
/^^(K) = —Kir/V in the friction coefficient for active particles and negative value of j3e- 
Because there is also normal positive friction obliged to the surrounding medium Pq the 
total friction is the sum /3s = /5o + /^^(V^)- For the simplest case is V-independent. For 
more complicated models it can be a slow function of y, as it was shown in [3]. We will 
consider here only this simplest model. The stationary distribution function for the active 
particles is: 



where = Dq -\- is the total diffusion, which is a constant if is V-independent. 
The coefficients in Eq. (49) are found on the basis of a kinetic consideration [3]. The 
characteristic parameters determined above in the case under consideration can be taken 
= /5o, Vq = Kgr/Po, Dc = Dy, and d(y) = 1. The physical reason for the maximum of 
the distribution function at V = K^/ISq is evident. It is due to the balance of the normal 
friction and self-accelerating forces. For self-motion of grains (cells) for example in a gas of 
hard spheres we can use for /3o and Dq the expressions (15, 17). A similar ^-dependence of 
the distribution function was found experimentally and phenomenologically for cells [12,14]. 
Using Eq.(43, 50) we find (for s = 3): 



In many other cases, when the microscopic theory of active motion is not developed yet the 
velocity independent diffusion also seems the suitable approximation [15]. For the all these 
models Eq. (42) can be simplified: 



If the characteristic constant value for the function 13{P) is j3c and the characteristic mo- 
mentum is Po (e.g. for the active particles the value of momentum, where /3(Po) = 0) we 
can introduce the undimensional function (f){P/Po) by the equality /3{P) — l3c(t>{P/Po) and 
rewrite Eq. (42) in the form: 




(49) 




(50) 




(51) 



Dr 




0(A), 



(52) 



sM^D 
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where the parameter A = P^Pc/MD and O, = Y/N is the the ratio of the integrals: 



/pT poo r 

d^T T dr^ d^^f^iO exp A / 
Jo Jri I Jt 



du}(j){uj) 



(53) 



N^Jd^CfM (54) 



For the self-consistent scheme the distribution function itself, according to Eq. (7), is equal 
to: 



/°(0=Cexp 



-A 

'0 



[ dp(j){p) 
Jo 



(55) 



VI. DIFFUSION IN THE COORDINATE SPACE ON THE BASIS OF 

MASTER-TYPE EQUATION 

Consideration of diffusion in the coordinate space and relation between the diffusion in 
coordinate and velocity space have been done in Section 3 on the basis of the Fokker-Planck 
kinetic equation. At the same time this approach is not enough for consideration of some 
special types of space diffusion, when the PT function has a specific form, in particular, 
posses the slow decreasing tails in the coordinate space, when the Fourier-components for 
such PT functions are absent. There are many very important processes, where different 
types of the anomalous diffusion exist [16]. This anomalous diffusion can be related not 
only with the coordinate, but also, in general, with time-dependent PT function. Usually 
the problem of anomalous diffusion is considering on the basis of fractional differentiation. 
This basis provides the universal description for probability distribution of grain at large 
distances. At the same time to describe all the distances we suggest to use the approach of 
PT function, when the fractional differentiation is not present at all. Naturally the results 
for the power-type PT functions are the same. But in general case, the PT function approach 
gives general and simple description for all interesting cases with a wide class of short-range 
behavior and a long distance tails in the coordinate space. We also formulate the time- 
dependent equation for anomalous diffusion, which gives useful general description for many 
applications. These applications will be considered in details in the separate paper. 

Let us consider the diffusion in the coordinate space on the basis of master equation, 
which describes the balance of the grains incoming and outcoming the point r in the moment 
t. The structure of this equation is formally similar to the master equation in the momentum 
space Eq. (6). Naturally, for the coordinate space there is no the conservation law, similar 
to the momentum space: 



dt 



j dv' {W{v, r')f,(r', t) - W(r', r)/,(r, t)} . (56) 



The probability transition W{r, r') describes probability for a grain to transfer from the 
point r' to the point r per unit time. We can rewrite this equation in the coordinates 
p = r' — r and r as: 
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dt 



J dp {Wip, r + p)/,(r + p, t) - W{p, r)/,(r, t)} . (57) 



If we suggest that the characteristic displacements are small and expand Eq. (57) we arrive 
to the " Fokker- Planck" form of the equation for the density distribution fg{r,t) 



dfgiv,t) d 



dt drr, 



^„(r)/,(r,t) + J-(i?„^(r)/,(r,t)) 



(58) 



The coefficients and Bap describe the acting force and diffusion, respectively, can be 
written as the functionals of the probabihty function (PT) in the coordinate space W in the 
form: 

Aa{r) = J d'ppaW{p,r) (59) 

and 

^a/3(r) = ^ J d'pPaPpW{p,v). (60) 

For the isotropic case the the probabihty function depends from r the modulus of p. For 
homogeneous medium, when r-dependence of the PT is absent the coefficients — and 
the diffusion coefficient is constant Ba/s — SapB, where B is the integral 

This consideration cannot be applied to the specific situations, when the integral in 
Eq. (61) is infinite. In that case we have investigate the general transport equation (56). We 
will consider the problem for homogeneous and isotropic case, when PT function depends 
only from modulus \p\. By Fourier-transformation we arrive to the form of Eq. (56): 



^/p(k, t) 
dt 



J dV[exp(^kp)-l]iy(|p|)/,(k,t)=X(k)/,(k,t), (62) 



where ^(k) = X{k). Let us consider the simple form of PT function with a power depen- 
dence from distance W{p) — C/\p\°', where C is a constant. For one-dimensional case we 



find: 



X{k) = -4 r dpsm^ W{p) = -23-"C|fcr-i rdC^. (63) 



For the values 1 < a < 3 this function is finite and equal 



X(fe)^-2Cr(l-a)|fe|"-^cos ^-^ (64) 

where F is the Gamma-function. At the same time the integral (61) for such type of PT 
functions is infinite, because usual diffusion is absent. The considered procedure for the sim- 
plest cases of power dependence of PT function is equivalent to the equation with fractional 
space differentiation [17-19]: 
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^i^^CA^/'f,{x,t), (65) 

where A'*/^ is a fractional Laplacian, a linear operator, whose action on the function f{x) 
in Fourier space is described by A''/^/(a;) = —{k'^)^^'^f{k) = —\k\^f{k). In the considered 
above case /i = {a — 1), where < // < 2. At the same time for more general PT functions, 
which are (for arbitrary values p) are not proportional to the -ctth power of p, the method 
described above is also applicable, although the fractional derivative is not exists. 

For the cases of purely power dependence of PT the non-stationary solution for density 
distribution describes so-called super-diffusion (or Levy flights). The solution of Eq. (65) in 
Fourier space reads: 

fg{k,t)^eM-C\k\H), (66) 
which in coordinate space corresponds to a so-called symmetric Levy stable distribution: 

1 , 



Ob 

■,p,0 



{kty/^^ 



(67) 



{kty/^^' 

For general case as it follows from Eq. (62) 

fg{k,t)^CieMX{k)t], (68) 

with some constant Ci. 

Let us turn to the three-dimensional case for the PT function equal W — C/\p\'^ with a 
constant C > 0. Then the function X{k) is equal 



X{k) ^ in 1^ dpp' - l) H/(p) = 4^Cr-» jT" 



(69) 



This integral is finite for the values 3 < a < 5. If we return to the variable p = ot — s 
(s=3 for 3-dimensional case and equal to the power of k in the function X(A;)) we find the 
same limitation as in one-dimensional case < // < 2, which characterizes the fractional 
derivative. 

Naturally, consideration on the basis of PT function given above, permits to avoid the 
fractional differentiation method and to consider more general physical situations of the 
non-power probability transitions. Let us consider that for a simple example. Taking (for 
one-dimensional case) the PT function VF(p) in the form 

W(p)=cl^i^2hM. (70) 

with p > 0, we arrive (e.g. for one-dimensional case) to the function X{ky. 

X{k) = -23-"C|ifcr-^ r^^ {'^-^M-'^{^C/\k\T]}sin'C ^ _2'-^C\kr'Y{a/\k\P,a). 

Jo C° 

(71) 
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As easy to sec the function Y(a/\k\P, a) is finite for 1 < a < p + 3, because for the small 
values of distances for p > divergence is suppressed also for some powers a > 3. Simple 
calculation for a = 2 and p = 1 leads to the result, which cannot be found by the fractional 
differentiation method 

Y{a/\k\,2) = - - arctan(|A;|/(7) + In [l + k^/a^] . (72) 

The asymptotic behavior of the function X{k) for ^ oo is similar, as it follows from 
Eq. (72), to the case W{x) — c/p°'. Therefore the universal behavior of X(k) is provided by 
asymptotical properties of PT function for 1 < a < 3. 

For the normal Gaussian diffusion (see e.g. Eqs. (58), (61)) for arbitrary s the function 
is equal X{k) = —Bk'^. It leads to the Gaussian distribution in the coordinate space and 
to the usual time dependence of the square-mean displacement (r^) ~ Bt. For different 
types of anomalous diffusion this relation is violated and the dependence (r^(t)) has to 
be calculated on the basis of concrete anomalous distribution, which is determined by the 
respective equation for /^(r, t). As we already demonstrated and will show below the crucial 
function, which determines the process is the PT function W . 

We also would like to mention in the connection of the problem of generalized description 
of diffusion that the time dependence of the density distribution for the time dependent PT 
function can be very different from the classical diffusion. This kind of problems relates to 
the class of the stochastic transport, which is very popular in modern physics and describes so 
called subdiffusive behavior [20], e.g. photoconductivity in strongly disordered and glassy 
semiconductors or the resonance radiative transfer in a plasma [21], which posses many 
other applications. Here we only shortly formulate the problem, leaving it for the detail 
consideration in the next paper. For this purpose we formulate a more general transport 
equation for density distribution: 

/,(r, t) = /,(r, t = 0) + jyrj dr' {W{r, r', r, t - r)/,(r', r) - W{r', r,r,t - T)/,(r, r)} . 

(73) 

For the case of stationary PT function we return to Eq. (56). If there is no memory, but 
PT function is a function of current time r we arrive to the equation more general than 
Eq. (56), which describes the density evolution, with prescribed time dependence of PT. In 
particular in that case for slow space dependence of the function W we find the equation 
for diffusion similar to Eq. (58)-(60) with time-dependent coefficients Aa{r,t) and Bai3{r,t), 
which are calculated on the basis of the PT function W{p, r, t). 

If the system posses memory, then, in the simplest case, when the function W can be 
expanded in spirit of Fokker-Planck approximation in the coordinate space we arrive to the 
following form of Eq. (73) 



/,(r,t) = /,(r,t = 0)+ fdr-,^ 







d_ 



A„(r, t - T)fg{r, r) + — (5„^(r, t - r)fg{r, r)) 



(74) 



If the function W posses memory and depends only from the difference t — t, but 
cannot be expanded in the coordinate space as in the case (74) we can use the Laplace- 
transformation in time to find: 
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/^(r,^) = M^l^_^ + J dr'{W{r,r',z)f,{r',z)-W{r',r,z)f,{r,z)}. (75) 

To arrive to this equation we suggested that due to fast convergence the upper hmit t in 
the integral (irTy(r) exp(— 2;t) with Im z > can be changed on infinity and determine 
W{r', r, z) by the equahty: 

POO 

W{r',r,z)= dTW{r',r,T)exp{-zT) (76) 
Jo 

(The same approach can be used to simplify the Fokker-Planck-type equation (74)). For 
the space homogeneous case Eq. (75) can be Forier-transformed and rewritten in the k, z 
variables: 

where 

X{k,z) = j dV[exp(ikp) - l]W{\p\,z). (78) 

If the PT function is time- independent W^(|p|,z) = W^(|p|)/2; and /c,(k, t = 0) =const. we 
return, as easy to sec, to the case of abnormal diffusion considered above in Eqs. (62)- (72). 
For a generale multiplicative form of PT function Vr(|p|,r) = W^i(|p|)M/2(t) the function 
X{k,z) is the product: 

X(k, z) = X^{k)X2{z) = J d'p [exp(ikp) - 1] Wi{\p\) ■ ^ dTW2{r) exp(-zT). (79) 

For the power r-dependence W2 — C/t^ with 7 < 1 the integral X2{z) is equal 

^2(^) = ^r(l-7). (80) 

Z ' 

The distribution function (77) in k, z space in that case reads 

^'^^^'^=[z-z^T{l-^)X,{k)y ^^^^ 

The detail consequences and physical applications of more complicated dependence 
H^(|p|, z) from z will be considered in the separate paper. 



VII. GENERALIZED FRICTION 

Let us consider now the generalized expression for the friction force Fy^ on the basis of the 
Fokker- Planck equation for charged grains. As it was shown in [3] for scattering of the ion 
stream with velocity u by moving grain with momentum Pq the of PT function ^^(Po, y, q) 
satisfies the equality: 
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w,(Po,y,q) = 'w;,(Po-y,q), (82) 

where ^^(P, q) is the PT for scattering of an isotropic system of ions by moving grains and 
y = Mu. The generahzed friction force for one grain can be represented through the friction 
coefficient: 

F^,(Po - y) = -(Po - y) /3.(|Po - y|). (83) 

In particular, Eq. (83) describes both hmiting cases: the friction force itself (for u = 0) and 
the drag force itself (for Pq = 0), acting on the grain due to ion scattering with a non-zero 
driven velocity. For two hmiting cases of the friction force itself Fjos for the grain with 
momentum G and immobile ions and for the opposite case - ion drag itself D^^ with the 
ion velocity u = G/M there is natural relation F/os(G) = — D^g(G). This picture leads 
to the (for practical applications important) generalization for a few species of the lighter 
particles, which we consider below. The particular result for D^^ for the case | u 
and the Coulomb scattering cross-section can be written in the form: 

D^5 = 2AoMurMnA, (84) 

where Aq — ^^(mi/M)a'^niVTi, = ^JT~[mi^ and a is the radius of grain. For dusty 
plasmas usually the parameter F = e^ZgZ^j oT^ ^ 1, and the Landau logarithm In A is given 
in [22], [23]. 

For the collecting PT function Wc(P,y,q), based on the existing collecting collision 
integral [7], a representation similar to Eq. (82) is not valid [3]. Instead of Eq. (82) we have: 

Wc(P,y,q) = Wc(P - y,q + ?^u), (85) 

where u'c(P, q) is the PT for the collecting process (in the case under consideration for the 
ions) in the isotropic ion system (u = 0), e.g., the Maxwellian one. The collecting part of 
the ion friction force for r\ = j2v\-^ <C 1 and arbitrary F has a form [2], [3]: 

F/e - -M\(3{rj, F) = -2M\Ao 1 - F - ^(1 - 3F) . (86) 

For F > 1 this formula describes a negative collecting friction F/c parallel to V. At the 
same time the collecting ion drag force ^ (derived in [22] and rigorously justified later in 
for ^ = V? /2v\-^ -C 1 and arbitrary F can be written as: 



Dj^ = 8AoMu 



F 

^+2 



(87) 



The difference of the modulus of forces in Eqs. (86) and (87) (for the case u = V) is the 
reflection of the collecting PT property, expressed by Eq. (85) and the structure of the 
collision collecting integral in dusty plasmas, when all the ions which cross the grain surface 
are absorbed and therefore completely transfer their momentum to the grain, although 
the mass transfer is considered as negligible. The problem of the mass transfer has been 
considered in [?]- [?] where a mass kinetic variable for dusty plasmas was introduced. 
The explicit solutions obtained in these papers were non-stationary. For RF dusty plasmas 
we are interested to find stationary solutions for the grain distribution function. These 
solutions take into account the surface recombination of the ions that bombarded the grain 
and remove the created atoms to ambient plasmas. Such a kinetic model (very close to many 
real experiments in RF discharge with dust) is in the process of being developed. 
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VIII. NEGATIVE GENERALIZED FRICTION FOR THE CASE OF ION FLOW 



Let us consider now the important generalization of the above approach for the case of 
two species of light particles, which transfer momentum to the grains. We assume an ion 
stream with velocity u, which can be scattered and absorbed by the grains (with velocity V) 
and an atomic subsystem with an isotropic velocity distribution function, e.g., a Maxwelhan 
distribution. Then by using the method developed in [3] and assuming that for the consid- 
ering model the equality (82) is satisfied, we arrive at the total force acting on the grains as 
expressed through the total friction coefficient PtiV): 

A^{P,y)^PMV), (88) 



(89) 



Using the approximations for atomic and ion friction coefficients f^aiV) — Pa{0)- where 
Pa{0) = 8^^{ma/ M)a'^naVT^ (scattering of a point atom on the hard sphere grain) and 
(3i{V) ~ 2Aor^lnA. This approximation for (3i{V) is valid for dominant ion scattering 
drag [24] and permits to avoid the problem of possible negative friction related with the 
mechanism of ion absorption, discussed in [2]. Such a way we arrive, for the simplest case, 
when the vectors of P and u arc parallel, at the picture typical for negative effective total 
friction. The total friction coefficient 



1 - 



2VmanaVT, 



(90) 



becomes negative for small grain velocities. 



IX. CONCLUSIONS 

The paper devoted to consideration velocity-dependent friction and diffusion coefficients 
in the velocity space and diffusion, including anomalous one, in the coordinate space. We es- 
tablished the general relations, which determine these coefficients (including tensorial struc- 
ture of diffusion) and showed, that this relations exist not only for the Fokker-Planck equa- 
tion, which applicability is restricted, at least for the Boltzmann-type of collisions, by low 
range of grain velocity. On the basis of master equation with the prescribed FT function we 
showed validity of these general relations for arbitrary grain velocities. 

We also considered the application of the PT function method to the ion absorption 
collisions in dusty plasmas, where due to the specific structure of the absorption cross-section 
and PT function the negative friction appears for a domain of low grain velocity. Necessary 
to mention that this result is valid for the considered model [2], but for application of the 
results to the real dusty plasmas, where the additional processes of surface ion recombination 
and atom remove to the plasma are exist, the elaborated models have to be developed. 

The important question about relation between the diffusion coefficients in the velocity 
and coordinate spaces is solved for the Fokker-Planck equation. For open systems, or for the 
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systems with velocity dependent coefficients the relation with coordinate diffusion is essen- 
tially more complicated, than in the usual case of time independent coefficients, connected 
by the Einstein relation. 

The problems of anomalous diffusion are considered on the basis of PT for the master- 
type equation in the coordinate space. Consideration presented above permits to avoid 
the method of fractional space differentiation (coinciding with the results obtained by this 
method in the particular cases of power-type dependencies of PT) and to extend the results 
for a wide class of PT functions. The extension of this approach is used to formulate the 
general integral equation for distribution function in the coordinate space, applicable to 
description of very different types of time-dependent normal and anomalous diffusion. This 
results, as we trust, can be important for very different applications. 

In many cases the PT function in the velocity space depends from some other vectors, 
except grain velocity. It can be ion driven velocity, for example, in dusty plasmas, or "driver" 
vector of cells in biological systems [3], or some additional vector for colliding inelastic 
grains etc. For these cases we have consider the methods of PT function determination and 
calculation the generalized friction and diffusion. Some development of these approach is 
also presented in the paper in application to scattering and absorption of ions by grains in 
dusty plasmas. 

For the case of scattering of ions, which have some driven velocity, by grains, moving in 
the slightly ionized dusty plasmas with essential friction between grains and atomic compo- 
nent we show manifestation of negative friction for low grain velocities. This process can be 
easily realized for some parameters of dusty plasma, providing acceleration and heating of 
grains. 
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